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ABSTRACT 

We construct the graviton propagator on de Sitter background in exact de 
Bonder gauge. We prove that it must break de Sitter invariance, just like 
the propagator of the massless, minimally coupled scalar. Our explicit so- 
lutions for its two scalar structure functions preserve spatial homogeneity 
and isotropy so that the propagator can be used within the larger context of 
inflationary cosmology, however, it is simple to alter the residual symmetry. 
Because our gauge condition is de Sitter invariant (although no solution for 
the propagator can be) renormalization should be simpler using this propa- 
gator than one based on a noncovariant gauge. It remains to be seen how 
other computational steps compare. 
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1 Introduction 



There are bound to be interesting lessons when two plausible arguments 
lead to opposite conclusions. The conflict we have in mind concerns the 
graviton propagator on de Sitter background and the arguments about it 
derive, respectively, from cosmology and mathematical physics: 

• From the perspective of cosmology, the graviton propagator cannot 
be de Sitter invariant because the unique Fourier mode sum for an 
invariant propagator is infrared divergent as a consequence of the scale 
invariance of the tensor power spectrum [1]. 

• From the perspective of mathematical physics, the graviton propagator 
must be de Sitter invariant because explicit, de Sitter invariant solu- 
tions for the graviton propagator equation have been obtained when 
covariant gauge fixing terms are added to the action [2J. 

The cosmology argument turns out to be correct. Seeing this teaches: 

• There is a topological obstacle to adding covariant gauge fixing terms on 
any manifold, and for any gauge theory, which possesses a linearization 
instability [1]. 

• One can impose covariant gauges in which the field obeys exact condi- 
tions, but previous solutions employed analytic continuation techniques 
that incorrectly subtract off power law infrared divergences |3]. 

• Solutions exist to the propagator equation which do not correspond to 
propagators in the sense of being the expectation value of the time- 
ordered product of two fields in the presence of some state 

The first point occurs even for flat space electromagnetism on the manifold 
X R: the invariant equations' linearization instability requires the total 
charge to vanish, whereas the Feynman gauge equations can be solved for any 
charge. The second point is familiar to everyone who has encountered the 
automatic subtraction of dimensional regularization. And a trivial example of 
the third point comes from multiplying the entirely real, retarded propagator 
by a factor of i. 

^It has been conjectured that this shows up in mathematical physics as a violation of 
reflection positivity [S]. 
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These insights resolve a number of puzzles in the literature. For example, 
employing the Feynman gauge fixing term for scalar quantum electrodynam- 
ics on de Sitter |6] produces a one loop self-mass-squared which possesses 
on-shell singularities |7]. These singularities seem to be the quantum field 
theory representation of what one would expect classically from an AqJ^ in- 
teraction energy in view of the erroneous temporal growth of Aq in Feynman 
gauge. The simplest noncovariant gauge [S] fails to show on-shell singular- 
ities 0. Nor is there any problem using the de Sitter invariant, Lorentz 
gauge propagator [9l[T0]. The conclusion for de Sitter electromagnetism is 
therefore that one must avoid adding covariant gauge fixing terms, but no 
physical breaking of de Sitter invariance occurs. 

The situation for gravitons is different owing to infrared divergences. It 
has long been noted that certain discrete choices of the two covariant gauge 
fixing parameters result in infrared divergences if one insists on a de Sitter 
invariant solution [111 112] ■ These choices had been dismissed as unphysical, 
"singular gauges" which must simply be avoided [13]. However, we can now 
see that they are precisely the cases for which the order of the omnipresent 
infrared divergence in the formal, de Sitter invariant mode sum changes from 
power law to logarithmic [3]. The power law infrared divergences of other 
choices were automatically — but incorrectly — subtracted by analytic reg- 
ularization techniques to produce solutions of the propagator equations that 
are not true propagators. The correct procedure in all cases is to allow free 
gravitons to resolve their infrared problem by breaking de Sitter invariance. 

The purpose of this note is to construct the graviton propagator in an 
allowed covariant gauge, without employing analytic continuation techniques 
to subtract off infrared divergences. Our procedure is to express the prop- 
agator in terms of covariant projectors acting on scalar structure functions, 
without making any assumption about the eventual de Sitter invariance of 
the result. These structure functions obey completely de Sitter invariant 
equations, but they fail to possess de Sitter invariant solutions on account 
of infrared divergences. The procedure is so general that we implement it 
as well for a vector particle of general mass My, and check that it agrees 
with the known de Sitter invariant solutions [3] for My > —2{D — l)if^ in 
the transverse sector and My > in the longitudinal sector. When de Sitter 
breaking must occur we have chosen to give explicit solutions which preserve 
the symmetries of homogeneity and isotropy that are relevant to cosmology. 
However, our equations for the structure functions are invariant, so one can 
easily derive solutions which respect any of the allowed subgroups. 
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Our notation is laid out in section 2. Section 3 presents a general treat- 
ment for minimally coupled scalars of any mass Ms- In section 4 we solve 
for the propagator of a vector with general mass My, including longitudinal 
and transverse parts. Section 5 applies the same technique to solve for the 
graviton propagator in de Bonder gauge. Our results are summarized and 
discussed in section 6. 

Because this work represents a long and mostly technical exercise we 
have thought it right to briefly discuss the physical motivation. The point 
is to facilitate the study of quantum effects from the epoch of primordial 
inflation for which the de Sitter geometry provides an excellent paradigmj^ 
The source of these effects is particle production. The small amount of 
particle production which has long been known to occur in an expanding 
universe [16] becomes explosive during inflation for any particle which is 
both massless and not conformally invariant on the classical level |T7]. This 
includes massless, minimally coupled scalars and gravitons |T8]. Of course 
this phenomenon is the origin of the tensor [19] and scalar [20] perturbations 
which are such an exciting tree order prediction of inflation. Our motivation 
is getting at the fascinating loop effects which should also be present. 

There have been extensive studies of the quantum loop effects from in- 
flation producing massless, minimally coupled scalars. In three models there 
are complete, dimensionally regulated and fully renormalized results: 

• For a real scalar with a quartic self-interaction there are one and two 
loop results for the expectation value of the stress tensor [21] and the 
self- mass-squared [22] • These show that inflationary particle produc- 
tion pushes the scalar up its potential, which increases the vacuum 
energy and leads to a violation of the weak energy condition on cosmo- 
logical scales without any instability. 

• For a massless fermion which is Yukawa-coupled to a real scalar there 
are one loop results for the fermion self-energy [23], for the scalar self- 

^Just how good can be quantified using the deceleration parameter q{t) — — aa/d^, 
which measures minus the fractional cosmic acceleration. Its value for de Sitter is g = — 1, 
and the threshold between inflation and deceleration occurs at g = 0. If one assumes single 
scalar inflation then the measured result for the scalar amplitude |14j , and the bound on 
the tensor-to-scalar ratio [11], imply 95% confldence that q{t) < —0.986 when the largest 
observable perturbations experienced first horizon crossing |15j . Because this would have 
been near the end of inflation, when q(t) was growing, most of the inflationary epoch was 
likely even closer to de Sitter. 
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mass-squared [21], and for the effective potential There is also 

a two loop computation of the coincident vertex function [25] • These 
results show that the inflationary production of scalars endows super- 
horizon fermions with a mass, which decreases the vacuum energy with- 
out bound in such a way that the universe eventually undergoes a Big 
Rip singularity [25] . 

• For scalar quantum electrodynamics there are one loop results for the 
vacuum polarization [26] and the scalar self- mass-squared [TJ [10] , and 
there are two loop results for coincident scalars [TD] , for coincident field 
strength tensors [27], and for the expectation value of the stress tensor 
[27] . These results show that the inflationary production of charged 
scalars causes super-horizon photons develop a mass, while the scalar 
remains light and the vacuum energy decreases slightly [28] . 

Scalar effects are simpler than those from gravitons because there is no 
issue about general coordinate invariance. They are also generally stronger 
because they can avoid derivative interactions. However, scalar effects are 
correspondingly less universal and less reliable because they depend upon the 
existence of light, minimally coupled scalars at inflationary scales. In four 
models with gravitons there are complete, dimensionally regulated and fully 
renormalized results: 

• For pure quantum gravity the graviton 1-point function has been com- 
puted at one loop order [2^. This result shows that the effect of infla- 
tionary gravitons at one loop order is a slight increase in the cosmolog- 
ical constant. 

• For quantum gravity plus a massless fermion the fermion self-energy has 
been computed at one loop order [30]. This result shows that spin-spin 
interactions with inflationary gravitons drive the fermion field strength 
up by an amount that increases without bound [3T| . 

• For quantum gravity plus a massless, minimally coupled scalar there 
are one loop computations of the scalar self-mass-squared [32] and the 
graviton self-energy [33]. The scalar effective field equations reveal 
that the scalar kinetic energy redshifts too rapidly for there to be a 
significant interaction with inflationary gravitons [32] ■ The effects of 
inflationary scalars on dynamical gravitons, and on the force of gravity, 
are still under study [33]. 
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• The nonlinear sigma model has been exploited to better understand 
the derivative interactions of quantum gravity p^l; explicit two 
loop results have been obtained for the expectation value of the stress 
tensor [35] . 

There are also a variety of other, sometimes less complete results, including: 

• In pure quantum gravity, a very early, approximate computation of 
the one loop graviton 1-point function was made [SB], as well as a 
later evaluation using adiabatic regularization [37]. Momentum cutoff 
computations of the one loop graviton self-energy [38] and the two loop 
graviton 1-point function [39j have also been done. 

• For gravity plus a scalar the one loop scalar contribution to the noncoin- 
cident (and hence unregulated) graviton self-energy has been computed 

m- 

• In scalar-driven inflation there have been computations of the one loop 
back-reaction [4lj, culminating in the realization that a physical mea- 
sure of the expansion rate shows no significant effect at one loop order 
[12] . (This issue is still open at two loop order, and as well in pure grav- 
ity [13].) There has also been a vast amount of work on loop corrections 
to the scalar power spectrum [U], including corresponding work on how 
to correctly quantify effects [45], and a powerful theorem by Weinberg 
which limits the maximum rate at which corrections can display secular 
growth [inj. Recently there has been renewed attention to the problem 
of untangling infrared effects from ultraviolet divergences [17] . 

• In gravity plus generic matter much interest has been devoted to the 
recent proposal by Polyakov [48J (following numerous antecedents [HI 
[36| [TH |39| HI]) that runaway particle production may destabilize de 
Sitter space [50] . 

2 Notation 

In the coming sections we shall study the de Sitter background propagators 
of three kinds of fields: minimally coupled scalars with arbitrary mass Ms, 
vectors with arbitrary mass My, and gravitons. The respective Lagrangians 
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are, 



--d^ipd^ipg 
1 



1 



-g - -M'sr^ 
1 



--d^A^d^A^g'^'^gP''^-- {D-1)H' + M^ A.A^g^^- 
1 



16nG 



R- {D-2){D-1)H' 



-9 



(1) 
(3) 



Here D is the dimension of spacetime, H is the Hubble constant (which gives 
the cosmological constant {D — 1)H^ = A) and G is Newton's constant. We 
make no assumption that the vector is transverse, although the form of its 
mass term in (jj]) obviously derives from partially integrating and commuting 
covariant derivatives in the Maxwell Lagrangian, and then adding a spuri- 
ous longitudinal kinetic term. The propagator of such a field appears in 
projection operators, even though the associated field cannot be dynamical. 

We define the graviton field as the perturbation of the full metric g^^u^x) 
about its background value 



(4) 



Once this definition has been made, there is no more point to distinguish- 
ing the background from the full metric so we drop the overbar and refer 
to the de Sitter background as simply g^^^x). Graviton indices are raised 
and lowered using this background field, for example, /i^^ = g^^hp^. Covari- 
ant derivative operators and other geometrical quantities are similarly 
constructed with respect to the background. Of special importance is the 
Lichnerowicz operator which, when simplified using the de Sitter result for 
the curvature Rp^/.^, = H^ig^^pgua - 9pa9up), takes the form. 



1 r 



+ - 



n + {D-i) 



rl 



gf^'^gP-^gPiPg-W 



H'. (5) 



Here and henceforth parenthesized indices are symmetrized, and □ = g^^^D^Dy 
is the covariant d'Alembertian operator. With the help of (|5]) we can express 
the free part of the gravitational Lagrangian 02]) in a convenient form. 



G 



SttG 



Surface Term] - -/i^^D'^''^"/ip,. 



^ + 0{h'). (6) 
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Much of our work will be valid in any coordinate realization of de Sitter 
space. However, when breaking de Sitter is necessary we shall always do so 
on the D-dimensional open conformal submanifold in which de Sitter can 
be imagined as a special case of the larger class of homogeneous, isotropic 
and spatially flat geometries relevant to cosmology. A spacetime point x'^ = 
X*) takes values in the ranges, 

- oo < < and - oo < x' < +00 for i = 1, . . . , (D-l) . (7) 



In these coordinates the invariant element is. 



ds'^ = g^ydx^dx'' 



— {dx^Y + dx-dx = a^^rj^ydx^dx^ 



(8) 



where rj^i, is the Lorentz metric and Qx = —1/Hx^ is the scale factor. 

Although important de Sitter breaking occurs, it turns out that the vast 
majority of our propagator is de Sitter invariant. This suggests to express it 
in terms of the de Sitter invariant length function y{x; z), 



y{x; z) = axa^H"^ 



x — z 



-{\x'- 



z \—ie 



(9) 



Except for the factor of ie (whose purpose is to enforce Feynman bound- 
ary conditions) the function y{x; z) is closely related to the invariant length 
i(x; z) from x'^ to z^^, 



y{x; z) 



4 sin' 



'{^-Et{x-z, 



(10) 



With this de Sitter invariant quantity y[x] z\ we can form a convenient 
basis of de Sitter invariant bi-tensors. Note that because y{x\z) is de Sitter 
invariant, so too are covariant derivatives of it. With the metrics g^uix) and 
gfiu{z), the first three derivatives of y{x; z) furnish a convenient basis of de 
Sitter invariant bi-tensors [7], 



dyjx; z) 

^x^^ 
dyjx; z) 

dz" 
d'^y{x; z) 
Qxi^dz^ 



\.x. 



H 1 



Hax(y6l + 2a,HA 
Ha,{y5l-2axHAx, 
H^a^a, {y6y, + 2a,HAxX - 2aJlHAx, - 2r], 



(11) 
(12) 

(13) 
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Here and subsequently we define Ax^ = rj^y{x—zY . Acting covariant deriva- 
tives generates more basis tensors, for example [7j, 

» - ^■^(2-.)«..(x) . (14) 

= HH2-y)a,^(.) . (15) 

The contraction of any pair of the basis tensors also produces more basis 
tensors [7], 

g..,,,^^ = «VW-^?=|^|^, (19) 

^ ' dxi'dzP dx^dz"" ' Ozp dz" ^ ' 

= 4i7V(-)-^^|^|^- (20) 

^ ^ ^ ^x^'^zP dx^dz"" ^ Sx'^ ^ ^ 



3 Scalar Propagators 

Scalar propagator equations play an important role in our analysis because 
our strategy is to enforce the de Bonder gauge condition, without making 
assumptions about de Sitter invariance, using covariant derivative projectors 
acting on scalar structure functions. The graviton propagator equation will 
then be used to infer invariant equations for these scalar structure functions. 
The point of this section is to review and systematize previous work [H |3] 
about how to solve such equations. We begin giving a general scalar propa- 
gator equation and explaining why infrared divergences for M| < preclude 
a de Sitter invariant solution. We review the two fixes in the literature, and 
then give a simple approximate implementation for our favorite one. The 
section closes with some powerful results for integrating propagators. 

One can see from ([T]) that the propagator of a minimally coupled scalar 
with mass Ms obeys the equation, 

n-Ml 



iA(x: z) 



(21) 
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The plane wave mode function corresponding to Bunch-Davies vacuum is 



u,y{x^,k) = \l —Tz ax ^ Hl^\—kx^) where u 

V 4:H 



D-1^2 Ml 

^ . (22) 



2 ^ m 

The Fourier mode sum for the propagator on infinite space is, 



+e{z'^-x'^)u^{x°,k)u^{z^,k)\. (23) 



(27r 



The result is de Sitter invariant when the integral converges 
zAf{x;z) 



hd~2 r(i2_i+^)r( 



(47r: 



r(f; 



D-1 D-1 D y 



2 

g^'^r(f-i) [^4^f-i ^^1 , 1 , 2/N 

(4.)f jy ^^^(2+^'2-^^2-2u) 



z/: — : 1 

' 2 ' 4 



(24) 



2 ' 2 

r(^+z/)r(^-z/)r(i-f) d-i y 



r(i+z.)r(i-z.)r(f-i) 



' 2'4^ " 



(25) 



iTT 2 



r(f)r(i-f) - 



r(f + z/+n)r(f-z/+n) „-f +2 



r(3-f +n) (n+1)! ^A' 
^ + v+n)T{^-v+n) ^yy 



(26) 



r(f+n)n! 

The gamma function r(^^^ — u + n) on the final line of ( ]26ll diverges for, 



M| = -A^(D-1 + A^)if=^ . (27) 



Its origin can be understood by performing the angular integration in the 
naive mode sum ( 123|) and then changing to the dimensionless variable r = 
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x|0(x°-^°)/7«(-'^a;°)i/«(-'^-s°)* + (conjugate) |, (28) 



2^71 



^^^'^ ^^'^ ^) * + (conjugate) | . (29) 



In these and subsequent expressions we define Ax = — That the 
divergence at (!27|) is infrared can be seen from the small argument expansion 
of the Bessel function and from its relation to the Hankel function, 

Hl^\x) = ^Ii^ffll^{e— J.(.)-J_.(x)}. (31) 

The small r behavior of the integrand fl29|) derives from three factors, the 
first being r^~^. The second factor takes the form, 

(-y^^HAxrY' ' Jn^iV^.HAxr) = —jy^Y^Ciiny . (32) 

^ J- n=0 

And the final factor from the Hankel functions is, 

^..( /^,)^a,( /^,)- . f W^^l f c.(„)r- , (33) 

Hence the small r expansion of the integrand has the form. 



1 ^^..^_2. , 



1v CO 



r 



oo 



2r(.)r(2.) ^z.-2-2.^^^^^)^2„_ (34) 



vrir(^)r(i.+i) 



n=0 
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The naive mode sum fl23l) is infrared divergent for 

D-2-2u<-l M| < . (35) 

However, there will only be a logarithmic infrared divergence, either from the 
leading term in or from one of the series corrections at n = A^, if one 
has, 

D-2-2v + 2N = -l = -N{D-1 + N)H^ . (36) 

This is precisely the condition fl27|) for the formal, de Sitter invariant mode 
sum fl26l) to diverge. 

The infrared divergence we have just seen was first noted in 1977 for 
the special case of Ms = by Ford and Parker [SI]. The appearance of an 
infrared divergence signals that something is unphysical about the quantity 
being computed. The correct response to an infrared divergence is not to 
subtract it off, either explicitly or implicitly with the automatic subtraction 
of some analytic regularization technique. One must instead understand the 
physical problem which caused the divergence and then fix that problem. As 
we will see, the fix involves breaking de Sitter invariance, which was realized 
in 1982 for the special case of Ms = |55]. Allen and Folacci later gave a 
rigorous proof that de Sitter invariance must be broken [56] . 

The divergence fl35|) occurs because of the way the Bunch-Davies mode 
functions fl2^ depend upon k for small k. The unphysical thing about having 
Bunch-Davies vacuum for arbitrarily small k is that no experimentalist can 
causally enforce it (or any other condition) for super-horizon modes. This 
has led to two fixes: 

1. One can continue to work on the spatial manifold R^~^ but assume 
the initial state is released with its super-horizon modes in some less 
singular condition [57]; or 

2. One can work on the compact spatial manifold T^~^ with its coordinate 
radius chosen so the initial state has no super- horizon modes [58] . 

We will adopt the latter fix. This makes the mode sum discrete, but the 
integral approximation should be excellent, and gives a simple expression for 
the propagator which differs from (123!) only by an infrared cutoff at k = H. 
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From the preceding discussion we see that the infrared corrected propa- 
gator iA(x; z) is just with the lower hmit cutoff at r = l/^Uxttz, 

ffD-2 i-oo Jd-a {Ja^az HAxt) 

,A(x-z) = -^i^ / ^^"^ ^ 



2^7T— {\^a^,HAxT)^ 
x|e(a;°-^°)i7«(^r)/f«(^r)* + ^^(^°-x°)(conjugate)|. (37) 

Of course we can express the truncated integral as the full one minus an 
integral over just the infrared, 



oo roo 



1 



dT= I dT- I ^"^"dr ^ iA(x: z) = iAf(x: z) + A]^(x: z 



(38) 

In this case it does not matter if dimensional regularization is used to evaluate 
both %A^{x\ z) and A[,^(x; z) because the errors we make at the lower limits 
will cancel. 

A further simplification is that A}^[x\ z) only needs to include the infrared 
singular terms which grow as a^ciz increases. These terms come entirely from 
the J-y parts of the Hankel function and they are entirely real, 

. _^2r(.)r(2.) .^^^^^^_^J._^{^.HAxT) 



J-.U -r J_JJ-r .(39) 



(47r)f r(z/+i) io (iyS^T^i/Axr)^ 

r2(l-z/) 

The final result is [59l [3] , 



(4^)f r(i!fi)r(^+i) 



n — 1 r ~ 1 

(ax.a.)-w-^ f /g. , V r 21^—2™ r4m 



where the coefficients CMnm are. 



1\N r(D-l 



r — ^ y ^ 2 



Ti^ + N + n-u) 



m\n\{N-n-2m)\ T{i^ + N-v) 

r(i^) Tjl-u) Tjl-u) 

T{^ + N-2m) r(l-z/+n + 2m) r(l-z/+m) ' ^ ^ 
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Of course there is no point in extending the sum over to values N > u — 
{^-^) for which the exponent of a^a^ becomes negative. Those terms rapidly 
approach zero, and they can be dropped without affecting the propagator 
equation because they are separately annihilated by □ — M^. 

It might be worried that the approximations made in deriving the infrared 
correction do violence to delicate consistency relations in quantum field the- 
ory, but this is not the case. For the Ms = scalar renormalization has been 
successfully implemented at one and two loop orders [TJ [TOl [211 1221 |23l [251 
[26l [271 [28]. Because the physical graviton polarizations obey the same mode 
functions as massless, minimally coupled scalars [H], one can also test the 
integral approximation with the graviton propagator. There is no disrup- 
tion of powerful consistency checks such as the Ward identity at tree order 
[60] and one loop [39]. Nor is there any problem with the allowed one loop 
counterterms [291 [301 [32l [33] . 

It is worthwhile to summarize these results in the context of a consistent 
notation. Consider a general scalar whose mass obeys M^/H'^ = (■^^)^ — 6^. 
Its propagator iAb{x; z) obeys the equation, 

□. + (,=ff'-(^)Vl.A.(x;.) = ^^^. (42) 

V 2 ^ J yj-g 

We define the final result for iA,b{x; z) as the limit as v approaches h of two 
functions which we wish to consider for general index v. The first term in the 
sum is zA„ f X'^ clS defined by expression f[2B]) . The second term is /S^^[x] z), 
as defined by expression fHU]) . except that the sum over is cut off at the 
largest nonnegative integer for which N < b — (-^^), with Al^{x; z) defined 
as zero for b < (^). Hence our final result is, 

i Afe (x; = hm liAf {x; z) + A]^{x; z)] . (43) 

We shall make significant use of four special cases for which a separate 
notation has been introduced: 

(^) ^ zAs{x;z) = B{y), (44) 



2 

i^Y^) ^ iAAix;z) = A{y) + 5A{a,,a„y) , (45) 

(^) ^ ^Aw{x■,z) = W{y) + 6W{a,,a,,y) ,m 

bM = ^^JiD-l)iD + 7) ^ tAMix; z) = M{y) + 5M(a„ a„ y) . (47) 



>w 
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Although the 5- type propagator is de Sitter invariant, its A-tjpe, 14^- type 
and M-type cousins have de Sitter breaking parts, 



6A 
6W 

5M 



k\ia.{a^a^) 



\bM-bA 



D-1 



ln(a^a^)(y-2) 




(48) 
(49) 



bM-bA 



bM-bA-1 



X 



(y-2) 
46a 



+ 




The constants k and kM are, 
HD~2 y{D-1) 



k = 



(4vr: 



k 



M 



AbAibM-l) 
_ H^~^ T{bM)r{2bM) 



(50) 



(51) 



(47r)f r(6A)r(6M + : 

The main, de Sitter invariant parts of each propagator consist of a few, 
potentially ultraviolet divergent terms (at y = 0), plus an infinite series. For 
the M-type propagator there are no cancelations with the de Sitter breaking 
terms: just replace z/ everywhere by 6m in expression (|26ll to find M{y) = 
iAf^^{x] z). However, there are cancelations when this replacement is done 
for the A-type and W propagators, 

Biy) 



4xf-i 



My) 



Wiy) 




(52) 



(53) 



(n-f + 2)(n-f + l)(n + l)! n{n-l)T{n + - 



(54) 
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And the D-depdendent constants Ai, Wi and W2 are 
T{D-l) 



W2 



+ (55) 



V{D + l) \ D + l 

r(f+i)j 2D 

T{D + l) 



(56) 



(57) 



A problem we shall often encounter consists of integrating one propagator 
against another. The result can be represented as the solution of a modified 
propagator equation, 



The solution is easily seen to be [H |3] , 

1 



(58) 



iAfec(x; z) 



iAc{x; z)-iAh{x; z) 



iAcb{x; z) . 



(62-c2)i72 

For the special case that the indices h and c agree one gets a derivative, 

1 d 



(59) 



iAbbix] z) 



2bH^ db 



iAbix: z) 



(60) 



We can obviously continue the process ad infinitum. For example, con- 
sider the case where the source is not a propagator but rather a singly inte- 
grated propagator. 



D-U2 



iAhrJx: z) = iArdix; z 



(61) 



The solution can be written in a form which is manifestly symmetric under 
any interchange of the three indices a, h and c. 



iAbcd{x]z) 



iAbdjx; z)-iAbc{x; z) 

{(f-c^)iAb{x] z) + {}p'-(P)iAc{x] z) + {c^-h'^)iAd{x] z) 

(62-c2)(c2-rf2)((;2_52)^4 



(62) 
.(63) 
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The case in which two of the indices are the same gives, 



2cH'^ dc 

And equating all three indices produces 

1 d 



1 d iAcc{x;z)-iAbc{x;z) 
lAbcix; z) = - 



iAbbbix; z) 



2bH^ db 
1 



iAbJx: z] 



c=b 



863 i/4 



-^tAbix; z)-b[-^ytAbix; z) 



(64) 

(65) 
(66) 



4 Vector Propagators 

One can see from ([2]) that the vector propagator obeys the equation, 



n-{D-l)H^-M, 



igf,pS^{x-z) 



-9 



(67) 



Note that we do not assume transversality; indeed, the full vector propa- 
gator cannot be transverse because the right hand side of equation fl67|) is 
not transverse. The first part of this section describes how to decompose 
the full propagator into its transverse and longitudinal parts, without mak- 
ing any assumptions about its eventual de Sitter invariance. Our technique 
is to express these parts using projectors formed from covariant derivative 
operators, acting on scalar structure functions. In the second part we derive 
a scalar equation for the longitudinal structure function and solve it using 
the techniques of section 3. In the final part we carry out the same analysis 
for the transverse structure function. The techniques employed here are a 
paradigm for the work of the subsequent section on the graviton propagator. 



4.1 Transverse and Logitudinal Parts 



The full vector propagator can be written as the sum of a transverse part 
and a longitudinal part. 



[x: z 



z) + i 



A^ 



X] z) 



(68) 



In previous studies [6l Ej the vector propagator was expressed as a linear com- 
bination of de Sitter invariant basis tensors like those introduced at the end 
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of section 2. Then the coefficient functions were chosen to enforce transver- 
sahty (or longitudinahty). This method is not open to us because we cannot 
assume de Sitter invariance for general mass My- What we require instead 
is a covariant decomposition which entails no assumption about de Sitter 
invariance. 

The longitudinal part is easy, 



[x: z) = 



d d 
dxf^ dzP 



SAx] z) 



(69) 



This expression is longitudinal for any choice of the longitudinal structure 
function J^l{x'i z). After much consideration we decided to express the trans- 
verse part as, 



x; z] 



Vf{x) X Vf{z) X 



x; z) X 



Tlpxixw) St{x]z) . (70) 



These symbols require explanation. The differential operator V^^ is defined 
by writing the Maxwell field strength tensor as F 



a/3 



fj, fi fi 



(71) 



Note that acting Vf{x) x V^^ (z) on any 4-index, symmetric function of 
X and z produces something with the right properties to be a transverse 
propagator. Of course some choices for the 4-index function give simpler 
final results than others! The best selection seems to be taking two of the 
four indices to be more covariant derivatives and the other two to belong to 
the section 2 basis tensor f[T5]) which gives a Lorentz metric in the fiat space 
limit. This corresponds to form f lTUj) with the definitions, 



QaK{x;z) 
TZl3x{x;z) 



D D 



2H^ Dx" Dz^ ' 
1 d'^y{x; z) 

2H^ dx^dz^ ■ 



(72) 
(73) 



4.2 Solution for the Longitudinal Part 

To derive an equation for the longitudinal structure function we take the 
divergence of the full propagator equation fl67|) . substitute relations fl68|) . 
f l69|) and flTOj) . and then commute the derivative to the left. 



■{D-1)H^-M. 



V 
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n^-{D-l)H^-M^]^n,SL{x;z 



V 



i6 (x—z) 



Hence we conclude, 



n,-M^]n,SL{x;z) 



i6^{x — z) 



From relation fH2]) of section 3 this implies, 



U^Siix; z) = -iAJx; z) 



for 



D-U2 Ml 



V 



The final solution for Sl follows from relations (1581159 

1 



z) 



iA 4 fx; z) —iAhix: z] 



iAAhix: z) . 



•D-i- 



(74) 
(75) 

(76) 
(77) 



(78) 

(79) 
and the 



We remind the reader of special case A with index Ba 
explicit expansion for iA^ix; z) given by expressions (145|) . ( 148|) and (!53!l . 



4.3 Solution for the Transverse Part 



Substituting our explicit solution fl7^ for the longitudinal structure function 
into the full propagator equation (!67|) allows us to derive an equation for the 
transverse part that was previously guessed [9], 



n-{D-l)H'^ 



-Ml 



[x: z 



x — z] 



d d 
dxi^ dzP 



iA 



A{x; z 



10) 



The most effective technique for solving this equation is to reduce each side 
of the equation to the standard transverse form (1701) . We then read off a 
scalar equation for the transverse structure function St{x; z), which can be 
solved by the methods of section 3. 

It is best to begin by establishing some important properties of the 
quadratic differential operator, 



p ^ = pa/3 p) 



5fn - D^D, 



(81) 
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We shall always contract P^^ into some vector T^, so it is possible to commute 
the final covariant derivatives to reach the form, 



n-{D-i)H' 



(82) 



It is tedious but straightforward to show that the covariant d'Alembertian 



commutes with P,,^, 



P/CT^ . 



(83) 



Note also that P^^ is transverse on both left and right, 

Df^P^Tp = = PfDpT . (84) 

P^''^ must therefore be proportional to transverse projection operator. The 
proportionality factor can be found by squaring. Comparing relations 
and (15^ implies, 



n-{D-i)H' 



n-{D-i)H^ 



T, 



5) 



The relevance of P^'^ is that it gives the differential operators in front of 
the general transverse form (170]) . 



Vf{x)xV:\z)xQ^,{x-z) 



^P/(x) X p;{z) 



(86) 



Substituting ( I70|) into equation ( IHOj) . and making use or relations ( l86i) and 
( 183|) . implies. 



n-{D-l)H^-M. 



x; z) 



= Vfix)xP;\z) X Q^^ix; z)[n-{D-l)H^-Mmn^,ST] • (87) 

We need next to consider what the d'Alembertian gives when acting on the 
factors to the far right. 



^^^^^ 'D^d^^^''^''^''- 



Recalling the definition ( [73|) of TZi3x{x;z), and making use of relation ([T] 
from section 2, gives two identities we shall use in this section and the next. 



D 



-TZb\{x\ z] 



mzpx{x] z) 



1 t .dy 
2^"^(^fe' 
-H^n^xix^z) . 



(89) 
(90) 
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Substitute these in (1551) and pass the single derivative back outside to obtain 

dy dSr 





d 


r dy 




■dz^ 


d 


r dy 


dx^ 


■dz^ 



-St 



d^y 



St 



dx^dz- 



-St + -Rf 



St 



-S. 



+ n 



s 



T ■ 



(91) 
(92) 
(93) 



The first term on the right of equation fl93|) is longitudinal. In view of relation 
we therefore conclude, 



n-{D-l)H^-M. 



[x\ z] 



Vf{x)xV;:\z)xQ^,{x;z) 



n 



I3X 



n-{D~2)H^-Ml 



Si 



(94) 



It remains to reduce the right hand side of (IHOl) to the standard form (JTO 
we have adopted for transverse bi-tensors, 



p 



(95) 



^y—g dx^^ dzP 
K^ix) X 'Pfiz) X Q„,(x; z)\npx{x; z)Vi{x- z)] . (96) 



This is easily accomplished by acting PJ^{x) x Pf{z) on both forms. Acting 
this operator on fl95l) and making use of relation fl8^ gives, 



Pf{x)xPf{z)z 



p 



iS (x — z) 



-9 

iS^{x — z) 



-9 



(97) 
• (98) 



= -2H^Vf{x) X V:\z) X Q^^ix; z) [gp^- 
= -2H^Vf{x) X V:\z) X Q^^{x-z)[npx{x-z 

Acting instead on (p6!) and making use of relations (185!) and (!93|) gives, 
P,^{x)xPP{z)t[,Pp]{x-z) 

= Vf{x) X V:\z) X Q„«(x;^)[□-(D-l)i72l'[7^;3A^ll , (99) 

= Vf{x) X V:\z) X Q„«(x;^) 7^;3A 



□ -(L>-2)/7^ Vi 



. (100) 
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Comparing expressions (lUSi) and fllOOl) implies, 



2i6 (x — z) 



n-{D-2)H^ Vi{x-z) = -2H 



Relation f l42|) from section 3 — with the special case of h 



(101) 

to infer, 
(102) 



n-{D-2)H'^\Viix;z) = -2HHABix;z) . 

Now apply relations (15811601) to finally obtain the structure function for the 
transverse projection functional, 

Pi(x; z) = -2HHABBix; z) . (103) 

We have now reduced the transverse propagator equation to the form. 



Vf{x)xV^\z)xQ^^{x;z) 



n 



D-{D-2)H^-M^ 



V^^{x)xV^\z) X Q^^{x;z) 



n 



■/3A 



-2HHAbb 



The transverse structure function therefore obeys, 

n-{D-2)H'^-M^]ST = -2HHAbb{x; z) . 
Again making use of relations (15811601) . our solution for it is, 

2H^ " ' 



.(104) 



(105) 



Si 



M2 



iAbb + 



2H' 



iAb — iAc 



where c 



D-3x2 Ml 



(106) 



5 The Graviton Propagator 

The previous section provides a model for the analysis of this section, except 
that we immediately specialize to gravitons which obey de Bonder gauge. 



(107) 



The first task is to express the propagator of such a graviton in terms of 
covariant projectors acting on scalar structure functions. With just a small 
extension of our previous results we can then commute the differential oper- 
ator to act directly on the structure functions. The final step is identifying 
the de Bonder gauge projection functional. 
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5.1 Enforcing de Donder Gauge 

In de Donder gauge fll07p the propagator must obey the gauge condition on 
either coordinate and its associated index group, 



X t 



X t 











(108) 
(109) 



Just as was the case for the vector propagator with the analogous conditions 
of transversahty and longitudinahty, we seek here to enforce (1108111091) by 
acting covariant projectors on scalar structure functions. It turns out there 
are two ways to do it, corresponding to the spin zero and spin two parts of 
the graviton propagator. 



[x: z) 



al3 



per 



,x: z) + i 



A 



per 



\x: z 



(110) 



The spin zero part of the graviton propagator is almost as simple as the 
longitudinal part of the vector propagator. It is a linear combination of 
longitudinal and trace terms from each index group. 



per 



[x: z 



V^,{x) X Vp,{z) 



So{x]z) 



111) 



The projector Vuv is 



V. 



fii/ 



9pu 
D-2 



n+2{D-l)H'' 



;ii2) 



Unlike the spin zero part of the graviton propagator, the spin two part is 
both transverse and also traceless within each index group. Recall that we 
obtained the key projector for the transverse part of the photon propagator 
by writing the Maxwell field strength tensor as F°'^ = V^^A^. We similarly 
define the projector V^^'^^ by expanding the Weyl tensor in powers of the 
graviton field C"^^^ = V^^^^h^"" + 0{h^). The final result takes the form [33], 



pi/ — ^pu £)_2 



•\a5 



{D-l){D-2) 



_g^SgM 



V 



;ii3) 
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The various pieces of this are, 



1 



1 r 
2 



V 



(114) 

,(115) 
(116) 



Acting V^^'^^{x) X Vp^^'^{z) on any eight index, symmetric function of x and 
z would produce a transverse and traceless tensor but, as with the vector, 
it pays to select a simple form. The best choice seems to be taking half the 
indices in the form of more covariant derivative operators, and the other half 
from two factors of the mixed second derivative (fT3l) of the length function. 



x; z) 



V^!^^\x)xV;^''^{z)xQ^,xQ. 



■ye 



n,3x'^s^<S2{x;z 



■ (117) 



We remind the reader of the definitions ( !72ll73l) of Qanix; z) and lZi3\{x] z). 

We close this subsection by giving the propagator equation. Acting the 
Lichnerowicz operator ([5]) on the graviton field and making use of the de 
Bonder gauge condition (11071) gives. 



-D 



pa 



4^ 



(118) 



This means the propagator obeys a relation of the form. 



n.,-2H' 



pi/^pa 



[x; z) - -gp^[x} 



'-9p.{pga)u X 



n.^ + 2{D-3)H 
i6^(x — z) / 



^ A 



ix; z] 



-9 



+ Other Terms ,(119) 



where the "Other Terms" make the right hand side consistent with the gauge 
condition. However, the right hand side of flll9p cannot be symmetric under 
the interchange of x'^ and z'^ (and their associated indices) because the left 
hand side of the equation obeys de Bonder gauge on z^ but not on x^. It is 
better to subtract off a term proportional to the trace. 



n,-2H' 



pumper 



[X] z) - -gpu[x} 



9pu 



D-2 



'2H' 



X 



n^ + 2{D-3)H^ 
{^)[U + 2{D-1)H\ 



a A 
a'-^pa 



OA 

a'-^pa 



pumper 



(x; z) + H^g^^{x)i 



^ A 



[x: z 



{x\z) 

'x;z) , (120) 
). (121) 
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It is easily checked that fll2ip obeys the de Bonder gauge condition on both 
x'^ and z'^. Hence the right hand side of the equation is symmetric under 
interchange of and z'^ and can in fact be guessed [3], 



[x; z) + H'^gf,^{x)i 



^ A 



9p,u9pa 



D-2 



W^ 



+ - 



^7] 



■ (122) 



Here i[^A^] is the full vector propagator for the tachyonic mass My 
—2{D — which obeys the equation, 



n + (D-l)H' 



lx:z 



i9t,pS^{x-z) 



-9 



(123) 



Recall from section 4 that it has the form given by equations ( l68l) . ( 169|1 
and ( 170|) . From equations ( 1791) and f ll06p we see that the longitudinal and 
transverse structure functions are, 



Sl{x; z) 
St(x; z) 



-iAamIx; z 
1 



D-1 



iAsBix; z)+iABwix; z) 



(124) 
(125) 



5.2 The Spin Zero Part 

To derive an equation for the spin zero structure function we simply take the 
trace of the propagator equation f ll22p . Tracing on the left hand side and 
making use of relations f lll01lll2|) gives, 



n^+2{D-l)H' 



a'-^pcr 



(x: z 



Vp.{z) 



Sf)(x: z] 



■ (126) 



Tracing the right hand side of f ll22p and making use of relations f ll24p and 
dSlEOD implies. 



n,+2{D-l)H' 



^ A 



x; z) 
2 gp„i5^{x-z) 



D-2 



-9 



~2DIDIiAm{x-z) , (127) 



-2Vf,a{.z)iAM{x]z) 



128) 



24 



The equation for Sq{x; z) derives from comparing expressions (I126p and 



□ + 2(D-l)if2 u + DH 



Sq{x]z) = -2(^jj—^yAM{x;z) 



Its solution follows easily from relations (15811661) . 

2iAMM(x; z)-2iAMwix; z) 



So{x; z) 



{D-1)H^ 



(129) 



2{jj—^yAwMM{x;z). (130) 



5.3 The Spin Two Part 

This is the most complicated analysis we shall have to make and it is greatly 
facilitated by the analogy with what was done for the transverse part of the 
vector propagator in section 4.3. Here, as for that case, the first step is to 
derive an equation for the remaining (spin two) part of the propagator by 
subtracting off the part we already have. We then establish some identities 
for a differential projector which comprises the exterior operators of the spin 
two part f lll7p of the propagator. These properties allow us to pass the 
d'Alembertian in the propagator equation through to act on the spin two 
structure function S2{x; z). Squaring this operator also allows us to express 
the right hand side of the propagator equation in the same form (I117P with a 
known structure function. Comparing the two sides of the equation leads to 
a scalar differential equation which can be solved by the techniques of section 
3. 

We derive an equation for the pure spin two part of the propagator from 
the full equation f ll22p by substituting the spin zero structure function (I130p . 
with definitions (Ilimil2]) . Now move everything known to right hand side 
to reach the form, 



1 r 



n.,-2H^ 



fJ.i^'-^pa 



[x: z) = I 



t^i^-^ pa 



9p.{p9o-)v 



9pu9pa- 



i6 



x — z] 



D-2 



+ 



. 9 'D~] 
D^^D^^t[,A^]{x; z) 
A^]{x;z 



[x;z), (131) 
jj 2 

+ (— — 7)'Pp,u{x) X Vpa{z)iAwM{x;z) 

]{x;z) 



+ DpD. 



2 \+DlDli,^ 



.Wl 



+ D:DIiIA^]{x; 



z 



(132) 



It can easily be checked that the right hand side of (11320 is transverse and 
traceless on each index group. We will eventually reduce this transverse- 
traceless projector to standard form. 



M'^-' pa 



(x: z 



X] 



pa 



[z] X Q^^ X Q. 



-ye 



(133) 
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However, it is best to first concentrate on the left hand side of the propagator 
equation (11321) . 

In analogy with the transverse projector Pj^ defined in equation flHTj) . we 
define the transverse-traceless projector, 



(134) 



We shall always consider this acted on a second rank tensor F^s- From the 
expressions f lll3mi6p which define V^^^^ it is straightforward but tedious to 
reach the form. 



-D^D.D^D^F^H + 



D-l 



D^D,-g^,n D'^D^F^p-UF^ 



+D^D^F,p + {D + 2)nF^, 



2g^,FJ^-2DF^, .(135) 



(Note the multiplicative factor of -D — 3 which derives from the fact that the 
Weyl tensor vanishes for = 3.) It is easy to see from 0135^ that P^/'' is 
traceless on both the left and the right, 

g^'Vj'Fps = = Pj\gpsF) . (136) 

It is also transverse on any index, both on the left and the right. 



135 



=D 



•(p 



135 



135 



{DpFs)= =Pj\DsFp). 



(137) 
(138) 



These two properties are very important because the only terms in expression 
(11351) which don't involve either divergences or traces are. 



fjLU 



1 / D-3 
'2^D-2 



n-2H^' 



.(139) 
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Hence squaring P^/*^ gives 



1 /D-3 



a7 



2 

1 /L'-S 



□ -2if^ 



n-DH' 



Vj'Fps , (140) 



2 VD-2 



p /3<5 



U-DH' 



Fps ■ (141) 



We note in passing that the covariant d'Alembertian commutes with P^/'^, 
just as it did for the transverse projector P^f- 

Of course the relevance of the transverse-traceless projector P^,^^'^ is that 
two factors of it give the exterior operators of the spin two part of the prop- 
agator, 



fJ.i^'-^pa 



1 



[x: z 



Vj\x) X V.^>^'f>{z)\npx{x;z)n,^{x-z)S2{x-z) 



From the fact that the d'Alembertian commutes with P^/^ we see, 



(142) 



[x: z 



-^Pj^(x) X PjHz) X - 

/^JJ4 pu \ J pa \ J 2 



n^xTZs^Si .(143) 



The next step is to pass the differential operator through to the structure 
function, making use of identities fl89ti90p from section 4, 



npxns^nS2 + 2g^\x] 



I3\ 



Dx'^ 



DTZ, 



Dx" 



dS2 

dx^ 



+2^"^(x) ~ 



[□7^^A)7^^<^ + -Rpxi^-Rs^) S2 ,(144) 



(3X /v<5(/> 

D 



n+2H' 



So 



Dx^ 



dy 
dz^ 



'R'5(h<S2 



+ 



D 



Dx^ 



n 



I3X 



dy 

dz^ 



When the external operators are contracted into this the terms on the final 
line of fll45p all drop by virtue of either transversality or tracelessness. Hence 
we have, 



n-2H' 



A^ 



(146) 
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It is now time to reduce transverse-traceless projection functional f ll32p to 
standard form (11331) . Just as we did with the transverse projection functional 
of section 4, this is accomplished by acting P^^'^{x) x P^/"^(z) on both forms. 
When acting on expression fll32p tracelessness or transversality make all but 
the first term drop out, 



p^.r(x) X p 



per 



Z) I 



p2 



x; z) 



i6^{x — z) 



TZ TZ 



i6^{x — z) 



-9 



, (147) 
• (148) 



On the other hand, acting the same operator on (11331) . and making use of 
relations f HiTjl and f lT46|) . tells us, 

1 



Z] I 



p2 
Ml'-' per 



\x\z] 



P ( r 



X P 



pa. 



TZ IZ -I 



/L)-3x2 2 




to 


( ) c 


n-{D-2)H^ 


^2 



(149) 



Comparing f ll48p with f ll49p we infer an equation for the structure function 
of the transverse-traceless projection functional. 



n-{D-2)H^fv2{x;z) = 16H^[ 



D-2^2 i5^{x-z) 
* X 



D-3 



-g 



(150) 



The solution is easily constructed using relation (^2]) and successive applica- 
tions of dSlEo]), 



/ 4 \ 2 

^2(x;z) = (^) 



iAAA{,x] z) - 2iAAB{x; z) + iABB{,x] z) 



(151) 



The long-sought equation for the spin two structure function derives from 
the substitution in equation (I13ip of relations (I146p and (I15ip , 



-US2{x-z) = [——) 



iAAA{,x; z) - 2iAAB{x; z) + iAsBix; z] 



(152) 



2 ^D-?>- 
The solution can be found using relations f l6ni66|) from the end of section 3, 
32 



S2{x; Z) = 



iAAAA^x; z)-2iAAAB + iAABB{x] z] 



(153) 
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6 Discussion 



We have constructed the graviton propagator on de Sitter background in ex- 
act de Bonder gauge fll07p . Our result takes the form (11 101) of a spin zero 
part and a spin two part. Both parts are represented in terms of covariant 
differential projectors which automatically enforce the gauge condition, act- 
ing on scalar structure functions. Our form for the spin zero part is given by 
relations (lllip and (I112p . The spin two part (11171) has a complicated defini- 
tion involving relations (1113111 16p and (I72ti73p . By taking appropriate traces 
and commuting differential operators we eventually derive scalar equations 
(11291) and (11521) for the structure functions of the respective parts. These 
equations are then solved using the general scalar techniques explained and 
summarized in section 3. 

We emphasize that our forms for the spin zero and spin two parts of the 
propagator involve no assumption about de Sitter invariance, nor specializa- 
tion to any particular portion of the de Sitter manifold. The equations (11291) 
and (I152p we derive for the two structure functions are scalar equations, valid 
in any coordinate system and with no inherent assumption about de Sitter 
invariance. To emphasize this, we act extra derivatives so as to make the 
source on the right hand side proportional to a delta function in each case. 



It happens that neither the spin zero structure function (I130p nor its spin 
two counterpart (I153p is de Sitter invariant. For the spin zero case this is 
obvious from the presence of tachyonic mass terms in both of the differential 
operators on the left hand side of equation (I154p . The mass M| = —DH^ 
includes a logarithmic singularity which shows up even in analytic regular- 
ization techniques. For the spin two equation (I155P the squared operator has 
positive mass-squared M| = {D — 2)H^ and would not lead to breaking of de 
Sitter invariance were it alone. However, the cubed operator is the same as 
that for a massless, minimally coupled scalar — as might have been expected 
from Grishchuk's old result |18]. Allen and Folacci long ago proved that this 
has no de Sitter invariant solution [56] . 

Exact de Bonder gauge is interesting because de Sitter invariant construc- 
tions based on analytic continuation methods had previously dismissed it as 
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an infrared divergent special case [13]. In fact, all valid gauges show infrared 
divergences. The special thing about de Bonder gauge is that some of its 
infrared divergences are logarithmic so that they are not automatically (and 
incorrectly) subtracted by analytic continuation. In all cases the right way 
to resolve the infrared divergence is by breaking de Sitter invariance. 

We have gone to considerable lengths — in previous work [H |3] and again 
in section 3 — to elucidate precisely what goes wrong with previous con- 
structions [2] which seemed to produce de Sitter invariant results. However, 
it worth pointing out that the fact of de Sitter breaking was already obvi- 
ous to cosmologists from the scale invariance of the tensor power spectrum, 
which becomes exact in the de Sitter limit [S]. It was also obvious from the 
explicit form of a propagator constructed by mode sums on the open sub- 
manifold (for which there is no linearization instability) [HI [8] . On the open 
submanifold the \D[D + 1) elements of the de Sitter group break down into 
four parts: 

1. {D — 1) spatial translations; 

2. i(D — 2)(D — 1) spatial rotations; 

3. A single dilatation; and 

4. [D — 1) spatial special conformal transformations. 

The gauge condition only breaks the last of these, but the solution for the 
propagator additionally breaks dilatation invariance [6T| [8] . The physical de 
Sitter breaking of this propagator was demonstrated by Kleppe, who aug- 
mented a naive de Sitter transformation by the compensating gauge trans- 
formation needed to restore the gauge condition [62j. Had the propagator 
been physically invariant this technique would have revealed it. 

We should also comment on the apparent conflict of our result with the 
pro-invariance argument given by Marolf and Morrison [63], based on work 
by Higuchi ^64j. They dealt with free dynamical gravitons in a noncovariant 
gauge on the full de Sitter manifold and they were able to construct the 
complete panoply of mode solutions and inner products. This should imply 
a vacuum which is physically de Sitter invariant — that is, invariant once 
the compensating gauge transformation is included. We know of no problem 
with this work but it should be noted that the propagator one gets using only 
dynamical gravitons (that is, the spatial, transverse-traceless polarizations) 
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is not complete. It is like the purely spatial and transverse photon propagator 
of flat space electrodynamics in Coulomb gauge. To fully describe electro- 
magnetic interactions also requires the instantaneous Coulomb interaction. 
Both of these are part of the same propagator in a covariant gauge such as 
the one we employ here. 

The constrained, spin zero part of our propagator — which is missing 
from the transverse-traceless part — provides the largest source of the de 
Sitter breaking we found. It is relatively simple to show that the de Sitter 
breaking terms in So{x; z) do not drop out when acted upon by the spin 
zero projector V^y{x) x Vpa{z). The spin two structure function contains less 
severe de Sitter breaking terms of the form, 



3 

breaking h = \ 



1 r i'^ 

'52(a;;2)Jdc sittc. = ^Sfc[ln(a^a^)J . (156) 



It is possible that these drop out from the spin two part of the propagator 
(11171) after all eight of the derivatives have been taken. In that case our work 
would be fully consistent with that of Marolf and Morrison. However, what 
we expect is that one of the infrared logarithms survives, which seems to be 
indicated by the scale invariance of the tensor power spectrum. 

The fact of de Sitter breaking in this system cannot be disputed, but 
there is wide freedom as to how one chooses to manifest that breaking. This 
freedom amounts to picking the initial state. We have chosen the explicit 
solutions of section 3 so as to preserve the symmetries of homogeneity and 
isotropy, which allow one to view de Sitter as a special case of a spatially flat, 
FRW geometry. This choice is known in the literature as the "i?(3) vacuum." 
Readers who prefer to preserve another subgroup can do so by starting from 
our scalar equations (I154p and (11551) . 

We wrote this paper to help resolve the long-standing controversy about 
de Sitter breaking for free gravitons, however, it has other applications. One 
of these is to test for gauge dependence in quantum gravitational loop correc- 
tions from primordial inflation. Of course gauge-fixed Green's functions will 
show such dependence, mingled with valid physical information. In flat space 
we would sift out the gauge dependence by forming the S-matrix. That ob- 
servable is not available in cosmology [BS] , and there is not yet any consensus 
for what replaces it. One technique is simply to carry out computations in 
different gauges. It may be that the leading infrared logarithm contributions 
(e.g., the one loop contribution to the fermion field strength from inflationary 
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gravitons [30]) are independent of the choice of gauge. Now we can test this 
conjecture using a completely different gauge from the one [HH E] employed 
in all previous computations. 

Our propagator should also make renormalization simpler because it pre- 
cludes the appearance of noninvariant counterterms. These complicated the 
analysis for previous computations [301 132] ■ It may also be that the gauge 
condition fll07p and the special properties of the differential projectors in our 
propagator make actual computations simpler. That turned out to be the 
case with the vector propagator in Lorentz gauge [9] for a variety of one and 
two loop computations [101 1271 [28] . 
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